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Abstract 

A  literature  review  of  recent  research  on  two-phase  flow  in  micro-channels  is  provided  in  this 
article.  Researches  on  the  micro-hydrodynamics  concerned  with  two-phase  gas-liquid  adiabatic  flow 
characteristics  in  both  circular  and  non-circular  micro-channels  are  discussed.  This  review  aims  to 
survey  and  identify  new  findings  obtained  from  this  attractive  area,  which  may  contribute  to 
optimum  design  and  process  control  of  high  performance  miniature  devices  comprising  extremely 
small  channels.  The  results  obtained  from  a  number  of  previous  studies  show  that  the  flow  behaviors 
in  the  micro-channels  deviate  significantly  from  those  in  ordinarily  sized  channels.  A  clear 
understanding  of  the  flow  behaviors  encountered  is  necessary  to  answer  the  fundamental  questions 
related  to  the  two-phase  flow  phenomena.  Similar  to  the  ordinarily  sized  channels,  it  is  expected  that 
the  flow  pattern,  pressure  drop  and  void  fraction  will  influence  the  two-phase  pressure  drop,  holdup, 
system  stability,  exchange  rates  of  momentum,  heat  and  mass  during  the  phase-change  heat  transfer 
processes.  It  is  obvious  that  publications  available  on  micro-channels  are  relatively  small  compared 
to  those  for  ordinarily  sized  channels  and  further  systematic  studies  of  two-phase  gas-liquid 
adiabatic  flow  in  micro-channels  are  required. 
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1.  Introduction 


In  the  past  years,  studies  on  two-phase  flow  and  heat  transfer  characteristics  in 
micro-channel  flow  passages  have  become  more  necessary  due  to  the  rapid  development  of 
micro-scale  devices  used  for  several  engineering  applications  including  medical  devices,  high 
heat-flux  compact  heat  exchangers,  and  cooling  systems  of  various  types  of  equipment  such 
as  high  performance  micro-electronics,  supercomputers,  high-powered  lasers,  and  so  on. 

In  addition  to  the  attractive  phenomena  caused  by  the  significant  role  of  the  surface 
effects,  this  emerging  field  enables  us  to  develop  powerful  miniature  devices  for  operation 
in  a  domain  which  seemed  to  be  unfeasible  in  the  past.  However,  it  is  obvious  that  a 
comprehensive  understanding  is  still  lacking  on  the  trends  and  parameters  dominating  the 
flow  behavior  in  mini-  and  micro -channels.  Thus,  fundamental  heat  transfer  and  flow 
characteristics  problems  encountered  in  the  development  and  processing  of  micro¬ 
electronic-mechanical  systems  (MEMS)  have  become  a  significant  challenge  for  the  design 
and  control  of  these  micro-systems.  According  to  the  requirement  of  new  solutions  for 
dissipating  heat  and  keeping  uniform  temperature  distributions  in  modern  devices,  a  clear 
understanding  of  the  major  mechanisms  affecting  micro-scale  heat  transfer  is  essential  for 
optimum  design  and  process  control  of  micro-systems. 

The  two-phase  flow  and  heat  transfer  characteristics  in  small  channels  such  as  micro¬ 
channels  and  mini  channels  are  likely  to  be  strongly  dependent  on  the  surface  tension 
effects  in  addition  to  viscosity  and  inertia  forces  resulting  in  significant  differences  in  the 
two-phase  flow  phenomena  between  ordinarily  sized  channels  and  small  channels.  Several 
investigators  have  proposed  the  criterion  for  a  small  channel  based  on  different 
dimensionless  parameters. 

Serizawa  et  al.  [1]  described  the  criterion  for  classification  of  micro-channel  proposed  by 
Suo  and  Griffith  [2]  as  follows: 


2 

D 


where  2  is  the  Laplace  constant  defined  by 


q(pl  -  pg ) 


Brauner  and  Moalem-Maron  [3]  proposed  the  different  criterion  as  follows: 
Eo^(2n)2, 


(3) 
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Nomenclature 


C  parameter  dependent  on  Reynolds  number 
C0  distribution  parameter 

D  channel  diameter,  m 

Z>B  bubble  diameter,  m 

Dq  equi-periphery  diameter,  m 

Dh  hydraulic  diameter,  m 

d P  pressure  drop,  Pa 

dz  unit  length,  m 

Eo  Etovos  number 

F  time  fraction 

/  friction  factor 

G  mass  flux,  kg/sm2 

g  gravitational  acceleration,  m/s“ 

j  superficial  velocity,  m/s 

KEavg  two-phase  mixture’s  average  kinetic  energy,  kg/ms“ 

L  length,  m 

Ti\  1  for  laminar  flow;  0.25  for  turbulent  flow 

n2  1  for  laminar  flow;  0.25  for  turbulent  flow 

Re  Reynolds  number 

s  channel  gap,  m 

U  average  velocity,  m/s 

V  rising  bubble  velocity  relative  to  the  liquid  phase,  m/s 
We v  Weber  number  for  vapor 

w  channel  width,  m 

X  parameter  defined  by  Eq.  (30) 

Xann  parameter  defined  by  Eq.  (37) 

x  flow  quality 

ot  void  fraction 

am  mean  void  fraction 

X  Lockhart-Martinelli  parameter 

t } b  two-phase  multiplier 

X  Laplace  constant 

/a  dynamic  viscosity,  kg/m  s 

v  kinematic  viscosity,  m  /s 

p  density,  kg/m 

o  surface  tension,  N/m 

Subscripts 

ann  annular  flow  regime 

B  bubble;  slug  bubble 

c  contraction 

G  gas  phase,  vapor  phase 

GS  gas  slug 
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I  gas-liquid  interface 

int  intermittent  flow  regime 

L  liquid  phase 

LO  all-liquid 

LS  liquid  slug 

liq  liquid  flow  regime 

s  smaller  channel 

T  two-phase  mixture 

UC  unit  cell 

VO  all-vapor 

vap  vapor  flow  regime 


where  Eo  is  the  Etovos  number  defined  by 


Eo 


9(Pl  ~  Pg)d2 

G 


From  Eqs.  (l)-(4),  o  is  surface  tension,  g  is  gravitational  acceleration,  D  is  channel 
diameter,  pL  and  pG  are  liquid  and  gas  densities,  respectively. 

Unfortunately,  it  is  likely  to  be  insufficient  to  identify  the  distinction  between  small 
channels  and  ordinarily  sized  channels  by  these  equations.  Some  arbitrary  channel 
classifications  based  on  the  hydraulic  diameter  Dh  have  been  proposed.  Recently, 
Mehendale  et  al.  [4]  employed  the  hydraulic  diameter  as  an  important  parameter  for 
defining  the  heat  exchangers  as  follows: 


•  Micro-heat  exchanger: 

•  Meso-heat  exchanger: 

•  Compact  heat  exchanger: 

•  Conventional  heat  exchanger: 


1  pm^Z>h^  100  pm 
100  pm  ^  Dh  ^  1  mm 
1  mm  ^  Dh  <  6  mm 
Dh  >  6  mm 


Based  on  the  evolution  of  small  flow  channels  used  in  engineering  applications,  however, 
Kandlikar  [5]  proposed  the  channel  classification  as  follows: 


•  Conventional  channels: 

•  Mini-channels: 

•  Micro-channels: 


Dh  >  3  mm 

200  pm  ^  Dh  ^  3  mm 

10  pm  ^  Dh  ^  200  pm 


Although  there  have  been  a  number  of  criteria  based  on  the  hydraulic  diameter,  the 
channel  classifications  available  cannot  relate  the  channel  diameter  to  the  fluid  flow 
mechanisms.  Further  investigations  should  be  conducted  to  meet  a  more  general  definition 
dealing  with  the  channel  classification. 

To  obtain  an  appropriate  design  and  process  control  of  the  various  engineering  systems 
employing  these  small  channels,  the  two-phase  flow  characteristics  as  well  as  heat  transfer 
mechanisms  in  various  small  channels,  especially  for  micro-channels,  should  be  clarified.  In 
the  past  decade,  there  has  been  a  relatively  small  amount  of  publications  available  for 
micro-channels  compared  to  those  for  ordinarily  sized  channels. 
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The  aim  of  this  paper  is  to  present  a  review  of  the  most  common  types  of  research  works 
done  on  two-phase  gas-liquid  adiabatic  flow  characteristics  in  micro-channels,  and  to 
propose  the  areas  available  for  further  investigations. 

2.  Two-phase  Gas-liquid  adiabatic  flow  characteristics  in  micro-channels 

Up  to  now,  there  have  been  publications  dealing  with  two-phase  gas-liquid  adiabatic 
flow  characteristics  in  micro-channels  with  various  geometrical  configurations.  However, 
the  two  most  common  types  of  such  configurations  employed  by  previous  researchers  are 
circular  and  non-circular  micro -channels. 


2.1.  Circular  micro-channels 


Adiabatic  two-phase  air-water  flow  characteristics,  including  the  two-phase  flow  pattern 
as  well  as  the  void  fraction  and  two-phase  frictional  pressure  drop,  in  circular  micro¬ 
channels  were  experimentally  studied  by  Triplett  et  al.  [6,7].  For  flow  visualization  results 
corresponding  to  their  first  paper,  the  flow  patterns  observed  were  bubbly,  slug,  churn, 
slug-annular  and  annular.  These  experimental  results  were  compared  with  other 
experimental  data  obtained  from  previous  investigators.  The  two-phase  flow  pattern 
transitions  generally  disagreed  with  the  lines  calculated  from  the  existing  models  and 
correlations. 

According  to  their  second  paper,  the  measured  values  of  void  fraction  and  two-phase 
pressure  drop  in  bubbly  and  slug  flow  regimes  were  found  to  be  best  predicted  by  the 
models  based  on  the  homogeneous  mixture  assumption.  In  the  annular  flow  regime, 
however,  the  experimental  void  fractions  and  pressure  drops  were  compared  with  the 
models  as  well  as  the  available  correlations,  with  poor  agreement. 

Kawahara  et  al.  [8]  used  a  circular  fused  silica  micro-channel  with  a  diameter  of  100  pm 
to  investigate  two-phase  flow  characteristics.  They  observed  many  different  flow  patterns 
simultaneously  in  the  circular  channel  under  any  given  flow  condition  tested.  Due  to  such 
simultaneous  appearances,  they  used  the  probability  of  appearance  of  each  flow  pattern 
along  with  the  time-averaged  void  fraction  obtained  from  each  flow  condition  to  define  the 
flow  regimes  as  shown  below. 


•  Slug-ring  flow  was  defined  when  the  probability  of  “thin  liquid  film”  was  greater  than 
that  of  “ring  liquid  film”  and  the  time-averaged  void  fraction  was  less  than  0.8. 

•  Ring-slug  flow  corresponded  to  when  the  probability  of  “ring  liquid  film”  was  larger 
than  that  of  “thin  liquid  film”  and  the  time-averaged  void  fraction  was  less  than  0.8. 

•  Semi-annular  flow  was  associated  with  flow  in  which  the  flow  alternated  mostly  between 
“liquid  alone”  and  “ring  liquid  film”.  The  time-averaged  void  fraction  was  greater  than 
0.8. 


•  Multiple  flow  included  “liquid  alone”,  “thin  liquid  film”,  “ring  liquid  film”,  “thick 
liquid  film”  and  “deformed  interface”,  and  the  time-averaged  void  fraction  was  less 
than  0.8. 

The  flow  patterns  defined  were  used  to  develop  a  flow  pattern  map  which  was 
subsequently  compared  to  other  flow  regime  maps  existing  for  two-phase  flow  of  air-water 
in  ~lmm  diameter  channels.  It  was  noted  that  bubbly  and  churn  flow  patterns  were 
absent  in  the  experimental  results  due  to  the  laminar  nature  of  liquid  flowing  through  the 
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micro-channel.  The  low  value  of  the  time-averaged  void  fraction,  even  at  high  gas  flow 
rates,  indicated  the  large  slip  ratios  and  weak  momentum  coupling  between  the  phases. 
The  friction  factor  in  single-phase  liquid  flow  and  two-phase  friction  multiplier  were  also 
analyzed.  The  single-phase  friction  factor  was  found  in  reasonable  agreement  with  the 
conventional  correlation.  The  two-phase  friction  multiplier  data  correlated  well  with  the 
Lockhart-Martinelli  correlation  based  on  a  separated  flow  assumption,  but  the  agreement 
between  the  data  and  homogeneous  flow  model  was  generally  poor. 

Chen  et  al.  [9]  experimentally  studied  the  characteristics  of  adiabatic  two-phase  flow  of 
nitrogen-water  in  glass  tubes  with  two  different  diameters,  1.0  and  1.5  mm.  Four  common 
flow  patterns  including  bubbly  flow,  slug  flow,  churn  flow  and  annular  flow  were  observed 
and  a  so-called  bubble-train  slug  flow  previously  reported  in  a  large  channel  under  micro¬ 
gravity  effects  [10]  was  also  observed  and  discussed  in  this  work.  The  random  manner  of 
the  number  of  bubbles  counted  in  the  bubble-train  slug  flow  pattern  was  attributed  to  the 
turbulent  nature.  Moreover,  increasing  liquid  superficial  velocity  tended  to  lower  the 
number  of  bubbles.  A  modified  drift  flux  model  was  found  to  be  an  appropriate  approach 
to  correlate  bubble  velocity  for  the  three  different  flow  patterns  covering  slug  flow,  bubble- 
train  slug  flow  and  churn  flow.  The  correlation  is  expressed  as 


UB  =  0.932/1",  (5) 

where  f/B  is  bubble  velocity  and  jT  is  superficial  velocity  of  the  two-phase  mixture. 

Furthermore,  the  experimental  void  fractions  were  well  predicted  by  the  use  of  this 
modified  drift  flux  model  in  conjunction  with  the  following  fundamental  equation: 


a  = 


Jg 

uB’ 


where  a  is  void  fraction  and  jG  is  gas  superficial  velocity. 

Serizawa  et  al.  [1]  investigated  the  visualization  of  the  two-phase  flow  pattern  in  circular 
micro-channels.  The  flowing  mixture  of  air  and  water  in  channels  of  20,  25  and  100  jam  in 
diameter  and  that  of  steam  and  water  in  a  channel  of  50  |im  in  diameter  were  conducted 
experimentally.  Two-phase  flow  patterns  obtained  from  both  air-water  and  steam-water 
flows  were  quite  similar  and  their  detailed  structures  were  described.  The  study  confirmed 
that  the  surface  wettability  had  a  significant  effect  on  the  two-phase  flow  patterns  in  very 
small  channels.  It  was  also  found  that  the  air-water  two-phase  flow  pattern  transitions 
generally  agreed  with  the  lines  predicted  by  the  Mandhane’s  correlation  and  the  cross- 
sectional  average  void  fraction  concerned  with  bubbly  and  slug  flows  followed  the  values 
obtained  from  the  Armand  correlation  [11]  for  an  air-water  system. 

Chung  and  Kawaji  [12]  performed  an  experiment  in  order  to  distinguish  two-phase  flow 
characteristics  in  micro-channels  from  those  in  mini-channels.  Four  different  circular 
diameters  ranging  from  50  to  526  jam  were  employed  to  examine  a  scaling  effect  on 
nitrogen-water  two-phase  flow.  The  results  including  the  flow  patterns,  void  fraction  and 
two-phase  pressure  drop  were  analyzed. 

For  larger  tubes  tested  (250  and  526  jam),  they  observed  flow  patterns  consistent  with 
those  obtained  from  mini-channels  having  diameters  of  around  1  mm.  The  corresponding 
flow  patterns  were  bubbly,  slug,  churn,  slug-annular  and  annular  flow.  The  time-averaged 
void  fraction  agreed  reasonably  with  the  Armand-type  correlation  available  for  mini¬ 
channels.  In  the  case  of  two-phase  pressure  drop,  it  was  found  that  at  a  given 
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Lockhart-Martinelli  parameter  the  two-phase  friction  multiplier  tended  to  change  with  the 
variation  of  mass  flux,  similar  to  that  appearing  in  mini-channels. 

For  relatively  small  channels  (49.5  and  100  pm),  only  slug  flow  was  observed  under  the 
experimental  conditions  tested.  The  void  fraction  and  pressure  drop  characteristics  were 
generally  opposite  to  those  of  larger  channels  previously  discussed.  The  relationship  between 
the  void  fraction  and  volumetric  quality  was  found  to  be  non-linear  and  hence,  the  Armand- 
type  correlation  became  unpredictable  for  these  channel  sizes.  In  addition,  the  frictional 
multiplier  was  independent  of  mass  flux.  Apparently,  diameters  between  100  and  250  pm 
seemed  to  be  in  the  range  causing  the  notable  change  in  two-phase  flow  mechanisms. 
Additionally,  they  modified  Garimella  et  al.’s  [13]  model  to  predict  the  two-phase  pressure 
drop  for  micro-channels  having  diameters  of  less  than  100  pm.  In  the  proposed  model,  they 
assumed  that  the  flow  consisting  of  unit  cells  (UC)  was  completely  intermittent  as  shown  in 
Fig.  1.  For  each  unit  cell,  there  were  two  regions,  including  liquid  single-phase  flow  region  (or 
liquid  slug),  which  had  no  entrained  gas  bubbles  and  the  region  occupied  by  the  two-phase 
flow  of  liquid  and  gas.  The  latter  region  was  formed  by  axially  symmetrical  gas  slug 
surrounded  by  an  annular  liquid  film  with  uniform  thickness  around  the  channel 
circumference.  Because  of  the  viscous  effect,  the  velocity  of  the  liquid  film  around  the  gas 
slug  was  assumed  to  be  significantly  low  compared  to  that  of  liquid  and  gas  slugs. 

In  the  single-phase  region,  the  pressure  gradient  in  the  liquid  slug  is  evaluated  from  the 
Darcy-Weisbach  equation: 


Fig.  1.  Schematic  diagram  of  unit  cell  for  model  of  Chung  and  Kawaji  [12]. 
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where  (dP/dz)LS  is  the  pressure  gradient  in  the  liquid  slug,  and  /LS  is  the  friction  factor  for 
the  liquid  slug.  The  friction  factor  for  fluid  flow  in  a  circular  channel  is  given  by 


64 

R?yS  ’ 


for  laminar  flow  (Re ls  ^2100), 


/ls  — 


0.3164 


for  turbulent  flow  (Reys  <  100, 000) 
in  a  smooth  channel. 


Eq.  (9)  is  known  as  Blasius  equation.  In  the  model,  it  was  noted  that,  for  the  transition 
from  laminar  to  turbulent  flow  (2100<ReLS<4000),  the  Blasius  formula  was  also 
employed  to  obtain  the  friction  factor.  The  Reynolds  number  of  the  liquid  slug,  Rey s,  is 
defined  as 


Reis  = 


PyUhsD 

Pl 


where  /ay  is  the  dynamic  viscosity  for  liquid  phase.  In  Eqs.  (7)  and  (10),  f/Ls  represents  the 
average  velocity  of  the  liquid  slug  and  is  expressed  by 


Uls=j^,  (11) 

1  —  a 

where  jy  is  the  superficial  liquid  velocity. 

In  the  two-phase  region,  the  average  velocity  of  the  gas  slug,  f/GS,  is  taken  to  be  equal  to 
that  of  the  gas  phase,  UG 

UG  s  =  UG  (12) 


and 


U  G  = 


.J_G 

Ot 


(13) 


The  pressure  gradient  in  the  two-phase  region,  (dP/dz)T,  is  given  by 


dPN 

dz 


— /gs 


Pg(^gs  —  U\  y 


2D 


B 


(14) 


where  /gs  represents  the  friction  factor  for  the  gas  slug  and  the  determination  of  /G s  is 
similar  to  that  of  /LS.  The  Reynolds  number  of  the  gas  slug,  RcGS,  is  given  by 


Rees  = 


Pg(Ug s  ~  Ui)Db 
Pg 


where  rg  is  the  dynamic  viscosity  for  gas  phase.  Ux  and  DB  in  Eqs.  (14)  and  (15)  represent 
the  average  velocity  of  the  gas-liquid  interface  and  bubble  diameter,  respectively.  U\  is 
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defined  by  Chung  and  Kawaji  [12]  as  follows: 


(dP/dz)T 

16jUL 


The  bubble  diameter  is  assumed  to  be 


Ob  =  0.900. 


The  unknown  variables  including  ReG s,  Oi?  and  (dL/dz)T  are  obtained  through 
Eqs.  (14)— (16).  Finally,  the  total  frictional  pressure  gradient,  (dL/dz)uc,  across  the  whole 
unit  cell  is  expressed  by 


dP\  Los/dO) 
d  zJjLuc  \d  z)lsLuc 


where  the  relative  lengths  of  the  gas  slug  (Lgs/Lxjc) 
approximated  by 


and  liquid  slug 


C^ls/^uc) 


are 


Luc  V^B  / 


Cls  _  |  _  CGs 
Luc  Cue 


Based  on  the  micro-channels  having  diameters  of  100  and  50  pm,  the  experimental  two- 
phase  pressure  drop  data  were  predicted  well  by  this  slug  flow  model. 

Coleman  and  Krause  [14]  carried  out  an  experiment  to  measure  R-134a  two-phase 
pressure  drops  due  to  contraction  in  circular  micro-channel  headers.  The  measured  two- 
phase  pressure  drops  were  considerably  high  compared  to  some  existing  correlations. 
However,  the  authors  suggested  that  the  use  of  model  developed  by  Schmidt  and  Friedel 
[15]  in  conjunction  with  the  momentum  density  model  was  in  reasonable  agreement  with 
the  measured  data. 

Two-phase  air-water  pressure  changes  in  the  circular  channel,  having  sudden  expansion 
and  contraction  sections,  were  experimentally  studied  by  Abdelall  et  al.  [16].  Two  stainless 
steel  tubes  were  connected  to  form  such  sections.  The  diameters  of  the  larger  tube  and 
smaller  tube  were  1.6  and  0.84  mm,  respectively.  The  two-phase  pressure  changes  were 
found  to  be  considerably  small  when  compared  to  those  predicted  by  the  homogeneous 
flow  model.  This  indicated  significant  velocity  slip  in  the  neighborhood  of  the  flow 
disturbance  resulting  from  the  area  change.  On  the  other  hand,  the  pressure  change 
data  were  well  predicted  by  a  simple  one-dimensional  (ID)  flow  theory  with  the  slip  ratio 
model  based  on  the  ideal  annular  flow  regime  with  minimum  entropy  generation.  They 
also  developed  an  empirical  correlation  of  two-phase  multiplier  for  the  flow  area 
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contraction,  </>lo,c: 

—  =  120(XjReLO,s)-a7,  (21) 

X 

where  Rey Q  s  is  the  Reynolds  number  with  the  total  flow  in  the  liquid  phase  in  smaller 
channel  and  the  Lockhart-Martinelli  parameter,  is  expressed  by 


where  x  is  quality,  py  and  /lig  represent  liquid  and  gas  viscosities,  respectively.  Liquid  and 
gas  densities  are  denoted  by  pL  and  pq,  respectively.  In  Eq.  (21),  the  subscript  LO 
corresponds  to  the  value  with  the  total  flow  in  the  liquid  phase  and  the  subscript  c  and  s 
correspond  to  contraction  and  the  smaller  channel,  respectively. 


2.2.  Non-circular  micro- channels 


Triplett  et  al.  [6,7]  also  conducted  the  same  systematic  experiments  as  discussed  in  the 
previous  section  to  study  gas-liquid  two-phase  flow  characteristics  in  semi-triangular  micro¬ 
channels.  The  results  obtained  were  similar  to  those  obtained  from  circular  channels. 

Xu  et  al.  [17]  investigated  an  adiabatic  co-current  vertical  two-phase  flow  of  air  and  water  in 
rectangular  channels  (12  mm  x  260  mm)  with  mini-/micro-gaps  of  0.3,  0.6  and  1.0  mm.  For 
channel  gaps  of  1.0  and  0.6  mm,  the  flow  patterns  observed  were  similar  to  those  found  in 
ordinarily  sized  channels  but  peculiar  flow  structures  which  were  different  from  typical  two- 
phase  flow  patterns  were  observed  in  the  channel  having  a  micro-gap  of  0.3  mm.  In  the  0.3  mm 
gap,  they  found  that  bubbly  flow  was  never  observed  and  new  flow  regimes  including  slug- 
droplet  flow  and  annular-droplet  flow  were  addressed.  For  the  new  flow  regimes,  it  was  observed 
that  liquid  droplets  were  always  adhered  to  the  wall  surface  and  such  droplets  were  blown  by  the 
moving  gas  phase.  They  also  presented  the  models  to  predict  the  flow  regime  transitions  in 
vertical  rectangular  channels.  The  models  for  bubbly  flow  to  slug  flow  transition  and  slug  flow  to 
churn  flow  transition  were  developed  based  on  Mishima  and  Ishii  [18]. 

The  transition  from  bubbly  flow  to  slug  flow  is  described  by 


(— !— 
\0.3Co 


0.35  ^/gDQ 


where  CQ  is  the  distribution  parameter  and  De  represents  the  equi-periphery  diameter  and 
is  given  by 


2  (s  +  w) 

? 

71 


where  s  and  w  are  the  channel  gap  and  width,  respectively.  In  Eq.  (23),  the  distribution 
parameter,  CQ,  for  rectangular  channels  is  given  by 


Co  =  1.35  -  0.35.  /—. 

V  P  l 

The  slug-churn  flow  transition  took  place  when 


m  • 


a>a 


(26) 
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The  void  fraction,  a,  is  given  by 


a  = 


Jg 


Co(/g  +7l)  +  ^ 


where  V  is  the  rising  bubble  velocity  relative  to  the  liquid  phase  and  is  expressed  by 

V  =  035dgE.  (28) 

The  slug  bubble  mean  void  fraction  am,  shown  in  Eq.  (26),  is  defined  as 

am  =  1  —  0.813V'75,  (29) 

where 


/ 

V  2d  ApL  B 


1)(/g  +,/l)  +  0-35  GdE 


Lb  and  Ap  are  the  slug  bubble  length  and  density  difference  between  the  two  phases, 
respectively. 

Based  on  a  void  fraction  of  0.75,  the  transition  of  the  annular  flow  took  place  when  gas 
and  liquid  superficial  velocities  satisfied  the  following  relationship: 


1-33/g  —  4/l 


\ 


(Pl  -  Po)9i>h  +  128CL(i^)'  ’pl/1  „|/2 

T“/g  ■ 


9.24Cg 


Ph 

VG 


-n  i 


Pg 


l  +  1-5fe) 


1/3 


where  CL  and  Cg  are  dependent  on  the  Reynolds  number. 

In  Eq.  (31),  vL  and  vG,  are  kinematic  liquid  and  gas  viscosities,  respectively.  nx  or  n2  is 
equal  to  unity  for  a  laminar  flow  and  nx  or  n2  is  equal  to  0.25  for  a  turbulent  flow. 
However,  the  models  developed  were  found  to  be  predictable  for  gaps  larger  than  0.6  mm. 
For  the  0.3  mm  micro-gap,  they  suggested  that  an  appropriate  theory  consistent  with 
peculiar  flow  regimes  observed  was  needed  to  develop  the  flow  regime  transition  models. 

Zhao  and  Bi  [19,20]  performed  experiments  to  study  two-phase  flow  patterns  in  vertical 
upward  air-water  flow  in  equilateral  triangular  channels  having  hydraulic  diameters  of 
2.886,  1.443  and  0.866  mm.  The  images  of  the  flow  patterns  were  taken  by  a  high-speed 
motion  analyzer  along  with  the  transient  pressure-drop  measurements.  The  flow  patterns, 
including  dispersed  bubbly  flow,  slug  flow,  churn  flow  and  annular  flow,  were  observed  in 
the  two  larger  triangular  channels.  These  observations  were  similar  to  those  generally 
encountered  in  the  ordinarily  sized  tubes  placed  in  vertical  direction.  However,  it  was 
interesting  to  note  that  a  new  flow  pattern,  namely  capillary  bubbly  flow,  was  identified  in 
the  smallest  channel  in  addition  to  the  typical  flow  patterns,  excluding  the  dispersed  bubbly 
flow.  Another  interesting  point  to  note  was  that  the  slug-churn  flow  and  churn-annular 
flow  transition  lines  shifted  to  the  right  as  the  hydraulic  diameter  decreased.  In  addition, 
the  previous  flow  regime  transition  models  being  available  for  upward  gas-liquid  two- 
phase  flow  in  the  ordinarily  sized  tubes  could  not  be  used  to  predict  the  flow  transitions. 

According  to  their  second  paper,  the  experimental  investigation  for  analyzing  gas 
velocity,  void  fraction,  and  pressure  drop  was  carried  out  in  the  same  experimental  setup. 
They  reported  that,  based  on  the  distribution  of  the  parameter  suggested  by  Ishii  [21]  and 
based  on  negligible  drift  velocity,  the  drift-flux  model  was  in  reasonable  agreement  with  the 
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experimental  data  for  gas  velocity  obtained  from  all  three  test  sections.  In  addition,  they 
obtained  a  correlation  for  predicting  the  void  fraction  and  single-phase  pressure  drop.  The 
void  fraction  correlation  developed  was  similar  to  that  suggested  by  Armand  and  Treschev 
[11].  It  was  found  that  the  Lockhart-Martinelli  correlation  was  in  good  agreement  with  the 
two-phase  pressure  drop  data  when  the  proposed  new  single-phase  friction  factor 
correlation  was  employed. 

Chung  et  al.  [22]  explored  nitrogen  gas-water  two-phase  flow  characteristics  in  a  square 
channel  having  a  hydraulic  diameter  of  96  pm  and  compared  the  results  with  those  for  a 
100  pm  circular  channel  employed  by  Kawahara  et  al.  [8].  Similar  results  including  pressure 
drop  and  void  fraction  data  were  obtained,  regardless  of  the  channel  geometry.  In  both 
channels,  the  Lockhart-Martinelli  correlation  appeared  to  do  rather  well  when  compared  with 
the  pressure  drop  data  while  the  void  fraction  data  were  found  to  be  in  poor  agreement  with 
the  homogeneous  flow  model  as  well  as  the  Armand-type  correlation.  They  also  discussed  the 
effect  of  channel  shape  on  the  transition  boundaries  for  flow  patterns  observed. 

A  flow  visualization  study  to  clarify  the  flow  patterns  of  vertical  upward  gas-liquid  two- 
phase  flow  in  rectangular  mini-channels  with  hydraulic  diameter  ranging  from  1.95  to 
5.58  mm  was  carried  out  by  Satitchaicharoen  and  Wongwises  [23].  Air-water,  air-20  wt% 
glycerol  solution,  and  air-40  wt%  glycerol  solution  were  used  as  working  fluids.  In  the 
experiments,  they  employed  various  rectangular  test  sections:  20  mm  x  2  mm, 
40  mm  x  1  mm,  40  mm  x  2  mm,  40  mm  x  3  mm  and  60  mm  x  2  mm  at  an  equal  length  of 
1  m.  The  flow  phenomena,  which  were  bubbly  flow,  cap-bubbly  flow,  slug  flow,  churn  flow 
and  annular  flow,  were  observed  and  recorded  by  a  high-speed  camera.  The  effects  of  gap 
size,  channel  width  and  liquid  viscosity  on  the  flow  pattern  transitions  were  discussed. 

English  and  Kandlikar  [24]  carried  out  experiments  to  explore  the  effect  of  surfactants 
on  two-phase  pressure  drop  of  air-water  adiabatic  flow  in  a  1  mm  square  channel.  The 
conditions  they  tested  were  based  on  the  operating  control  in  a  proton  exchange  membrane 
(PEM)  fuel  cell.  The  surfactant  Triton  DF-12  was  employed  to  reduce  the  surface 
tension  to  0.034  N/m.  The  application  of  the  surfactant  solutions  appeared  to  have  no 
effect  on  the  pressure  drop  due  to  many  possible  reasons,  such  as  characteristics  of  the  flow 
pattern  affecting  interaction  between  the  two  phases,  insufficient  surface  tension  for  an 
observable  change,  and  temperature  variation.  The  accuracy  of  existing  models  for  two- 
phase  pressure  drop  prediction  was  evaluated.  Based  on  the  Mishima  and  Hibiki  model 
[25],  English  and  Kandlikar  [24]  also  developed  a  new  model  applicable  to  low  mass  fluxes, 
high  mass  qualities  and  annular  flow  by  applying  the  Chisholm  value  for  laminar-laminar 
flow  of  5  instead  of  the  original  constant  of  21. 

Two-phase  flow  characteristics  including  flow  pattern  and  pressure  drop  of  ethanol-C02 
flowing  through  a  converging  or  diverging  rectangular  micro-channel  were  explored  by 
Hwang  et  al.  [26].  Different  from  previous  works  associated  with  the  flow  in  micro¬ 
channels,  they  observed  stratified  wavy  flows  at  the  entrance  of  the  converging  micro- 
channel.  Peculiar  flow  patterns  and  bubble  interactions,  which  were  due  to  the  presence  of 
pressure  and  acceleration  effects  in  the  variable  cross-section  channels,  were  also  discussed. 
For  the  pressure  drop,  they  found  that  the  channel  pressure  drop  in  the  diverging  micro- 
channel  was  substantially  small  in  comparison  to  that  in  the  converging  micro-channel. 
However,  the  two-phase  frictional  multiplier  plotted  against  the  Lockhart-Martinelli 
parameter  for  both  channels  was  independent  of  the  liquid  flow  rate.  This  result  was  in 
agreement  with  that  observed  in  nitrogen-water  two-phase  flow  in  circular  micro-channels 
investigated  by  Chung  and  Kawaji  [12]. 
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Due  to  different  flow  regimes  raised  in  the  same  micro-channel  under  a  given  flow 
condition,  Jassim  and  Newell  [27]  proposed  void  fraction  and  pressure  drop  models  based 
on  the  probabilistic  flow  regime  map.  This  modeling  technique  was  developed  in  such  a 
way  that  flow  visualization  data  were  incorporated  into  the  models  developed  for 
predicting  pressure  drop  as  well  as  the  void  fraction  for  air-water  at  20  °C  and  refrigerants 
including  R410A  and  R134a  at  10  °C  with  mass  flux  values  ranging  between  50  and  300  kg/ 
m  s.  A  test  section  with  6-port  micro-channels  of  1.54  mm  in  hydraulic  diameter  was  used 
to  obtain  time  fraction  data  for  different  flow  regimes  including  liquid,  intermittent, 
annular,  and  vapor  flows.  The  proposed  probabilistic  model  for  predicting  total  pressure 
drop  in  micro-channels  can  be  expressed  with  many  flow  regime  time  fraction  functions 
and  is  given  by 


+  Tint 


+  T 


vap 


+  F 


ann 


where  Tliq,  Fin t,  Tvap  and  Tann  represent  time  fraction  curve  fits  in  liquid,  intermittent, 
vapor  and  annular  flow  regimes,  respectively.  The  pressure  drop  for  single-phase  flow, 
including  that  for  liquid  flow  (dT/dz)iiq  and  that  for  vapor  flow  (dT/dz)vap,  is  determined 
from  the  conventional  relation: 


'dP\ 

dz  J 


=f—(— 

liq  or  vap  ^P  v^h. 


(33) 


where  G  is  mass  flux.  The  pressure  drop  model  for  intermittent  flow  assuming  a 
homogeneous  density,  pT,  and  employing  the  average  kinetic  energy,  KEavg,  is  given  by 


dTN 

dz 


int 


=  0.045  TKEavg, 
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The  annular  flow  pressure  drop  model  is  evaluated  as  follows: 


'd 


=  ^ 


vo1  dz  ) 


VO 


dz  /  ann 

9  • 

where  (/>y0  is  two-phase  multiplier  and  is  expressed  by 

4>\0  =  exp(-0.046Vann)  +  0.22[exp  (-0.0024Tann)  -  exp  (-7Vann)]: 


(34) 
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where  We v  in  Eq.  (37)  represents  the  Weber  number  for  vapor  and  is  defined  as 

[( xG)2/pG ] 


Wev  = 


<j/Dh 


(38) 


In  Eq.  (35),  the  vapor  only  pressure  drop  is  given  by 

1  G2 


d  P\ 

dz  J 


—  2/vo  n 
VO  Pg 


(39) 


The  probabilistic  void  fraction  model  proposed  is  also  expressed  with  many  flow  regime 
time  fraction  functions  as  shown  by 


^total  —  TfiiqOCiiq  F int^int  T  T^yapOCyap  T  E'ann0(ann. 


where  otnq  represents  the  void  fraction  for  liquid  only  flow  and  is  equal  to  0,  while  avap 
represents  the  vapor  only  flow  void  fraction  which  is  equal  to  1 .  The  intermittent  flow  void 
fraction,  aint,  and  annular  flow  void  fraction,  aann,  are  defined  as  follows: 


_  (0.833  +  0.167x)x(l/pG) 
a‘nt  (1  -x)(l/pL)  +x(l/pG)  ’ 


^ann 


1  + 


Pressure  drop  and  void  fraction  predictions  based  on  the  probabilistic  flow  regime  map 
were  found  to  be  in  accordance  with  the  measured  data  obtained  in  micro-channels. 
However,  care  should  be  taken  when  the  models  are  needed  to  be  applied  under  conditions 
and  geometries  other  than  those  in  which  the  probabilistic  flow  regime  maps  were 
established. 


3.  Conclusions 

As  discussed  above,  this  review  reveals  that  there  is  only  a  small  amount  of  publications 
available  for  the  two-phase  gas-liquid  adiabatic  flow  in  micro-channels.  The  existing 
literature  is  still  lacking  in  various  aspects,  including  geometry  of  mixing  chamber 
configuration,  fluid  properties,  channel  orientation  and  relevant  flow  models. 

It  is  obvious  that  the  status  of  current  researches  on  phase-change  heat  transfer  and  two- 
phase  flow  phenomena  in  micro-channels  is  still  in  infant  stage  and  a  great  deal  of 
systematic  investigations  will  become  substantially  more  necessary  to  meet  general 
conclusions  needed  for  the  appropriate  design  and  process  control  of  several  engineering 
applications. 
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